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ABSTRACT 

Context. We present an alternate formula for calculating gravitational time delay. 

Aims. We use this formula to reduce geodetic Very Long Baseline Interferometry (VLBI) data, taking into account gravitational effects 
within the solar system, and to test general relativity. 

Methods. The alternate formula was obtained by expanding the conventional formula in a Taylor series. We show that the gravitational 
delay can be split into several terms including a term due to the coordinate transformation and terms that are explicitly linked to the 
light deflection angle. 

Results. Our formula is compared numerically with the conventional formula, and difference in arrival times within 1 ps are found at 
1° from the Sun for a full range of baseline lengths. 

Conclusions. We conclude that the standard reduction of geodetic VLBI data for the effects of general relativity is equivalent to 
displacing the reference radio sources from their original catalogue positions in accordance with the classical light deflection formula 
across the whole sky. 

Key words, astrometry - proper motions - reference systems - methods: analytical 


1. Introduction 

With geodetic Very Long Baseline Interferometry (VLBI) pre¬ 
cise group delays - the difference in arrival times of radio waves 
at two radio telescopes ( [Schuh & Behrend|2Q12| ) - can be mea¬ 
sured, and very accurate radio source positions can be obtained 
(~ 40 microarcsec l |Ma et al.||2009] )). There are numerous as¬ 
trometric and geodetic effects that must be accounted for when 
reducing high-precision VLBI data. One of these effects is grav¬ 
itational time delay, which is caused by the gravitatio nal fields 
of massive objects along the observer’s line-of-sight ( [Shapiro| 
1964| \\961) and is responsible for the measured group delays. 
These massive objects also bend the path of light from distant 
objects, causing an apparent shift of their position in the sky, as 
predicted by general relativity ( |Einstein|1916| ). Within the solar 
system both the Sun and Jupiter are massive enough to produce 
these effects. The conventional equation for calculating gravita¬ 
tional delay is formulated in terms of the positions of the radio 
telescopes within the barycentric reference frame of the solar 
system dKopeikin|l990} [Eubanks et al.|1991t[Soffel et al.|1991 ; 

|Klioner|1991| ), rather than the baseline length between the radio 

telescopes. Moreover, the total VLBI delay formula includes a 
coordinate term to transform from the geocentric to the barycen¬ 
tric reference frame. 

We propose an alternate gravitational delay formula using a 
Taylor series expansion. In Sect. we show how the conven¬ 
tional formula can be split into a sum of several terms. All the 
terms are expressed as a ratio of the baseline length to the Earth 
barycentric distance. One of the terms links the gravitational de¬ 
lay and the well-known formula for the light deflection angle at 


an arbitrary elongation from the Sun. Another term cancels the 
coordinate term explicitly presented at the equation of the total 
VLBI delay. It is more convenient to formulate the effect of gen¬ 
eral relativity so that it is free from this complicating coordinate 
term. 

In Sect.j^a numerical comparison between the new and con¬ 
ventional formulae is presented for the cases of the gravitational 
field of the Sun and Jupiter. |Kaplan| ( |1998| ) compared the astro- 
metric difference in general relativity effects between the group 
VLBI delay and optical observations, where the observables are 
measured in angular quantities. It was claimed that the differ¬ 
ence is less than 1 microarcsec using some numerical algorithms. 
Our analytical approach shows that the two approaches are abso¬ 
lutely equivalent for all elongation angles (the difference is equal 
to zero), except for the small area around the gravitational body, 
which is a special case. 

In Sect. we show that the formulae for the post-post- 
Newtonian effect can be developed using the same approach as 
for the standard monopole light deflection. 

In Sect. we discuss the application of the new formula to 
estimate the proper motion of an extragalactic radio source in¬ 
duced by the solar gravitational held. 


2. Developing the conventional formula for 
gravitational delay 

In addition to the three classical tests of general relativity (GR), 
a fourth test - the delay of a signal propagating in the solar grav¬ 
itational held - was proposed by |Shapiro| ( |1964| ) and is known 
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as the Shapiro delay. The positions of astronomical instruments 
on Earth are referenced to the solar system barycentre, and the 
VLB I delay is equal to the terrestrial time (TT) coordinate time 
interval between two eve nts of the signal a rrival at the first and 
second radio telescopes ( [Petit et J^|201Q| ). The difference be¬ 
tween the two Shapiro delays as measured with two radio tele¬ 
scopes at the barycentric distances ri and r 2 from a body of mass 
M gives a gravitational delay, Tgrav, which must be taken into ac¬ 
count during the standard reduction of the high-precision geode¬ 
tic VLBI data ( [Petit et al.|2010| ), 


'^grav ~ 


iy + l)GM |ri| + (s ■ rQ 
” |r2l + (s-r2)’ 


( 1 ) 



Fig. 1. Angles (p, if/, 6, and A, originated from the position of gravita¬ 
tional mass (B), quasar (Q), and baseline vector (b). If the Sun plays 
the role of the gravitational mass, then the point B in Fig. 1 is also the 
position of the solar system barycentre. 


where G is the gravitational constant, s is the unit vector in the 
direction of observed radio source, c is the speed of light, and 
y is the pa rameter of t he parametrised post-Newtonian (PPN) 
formalism ( [Will[[1993] ), equal to unity in GR. The classical 
Shapiro delay formula includes the positions of the remote body 
and the receiver, but the differential delay Q only depends 
on the positions of the two receivers in the adopted reference 
system. Therefore, the formula 0 is valid for all objects 
regardless of their distance unless the curvature of the observed 
light front cannot be ignored. All extragalatic radio sources 
meet this condition. 


We introduce the baseline vector, b, as the difference be¬ 
tween the two barycentre radius-vectors of two antennas b = 
r 2 - ri. Then formula 0 may be re-written as follows (we recall 
that 7 = 1): 


^grav ~ 


2GM (r 2 ^ - 2(r2 • b) + b^)5 + (r 2 • s) - (b • s) 
c3 ” |r 2 | + (s • rj) 


( 2 ) 


and changing |r 2 | to r 2 and |b| to b for simplicity, one can obtain 
from Eq. 0 


2GM 


‘ grav 


-b cos if/-b cos (p+ ^ (1-cos^ if/) 
r2(l-cos 6) 


1 /-Z?cos0^-Z?cos^\^ 

2 V r2(l-cos6') ) 


( 6 ) 


Angle 6 is given here with respect to the second station, that is, 
6 = O 2 for the sake of simplicity. The difference between 6i and 
O 2 is discussed in Sect.j^ Similar expressions (without te rms that 
include have been obtained for geodetic VLBI (e.g. 


stein et aI.[I983t[HeIlings 1986]) and for SIM interferome 
sion model delays ( [Turysliev|2()09[ ). 

Lor the general case of no coplanar vectors (Lig.[T]) 

cos i// = - cos ip cos ^ - sin ^ sin Q cos A, 


Linkel- 


:ry mis- 


(7) 


Lor typical baselines, |b| « |r 2 |, and we can use the Taylor series 
expansion (I -l- I -l- | - y for v << I to modify (^ as 

follows: 


and formula ([^ for arbitrary angles p and 6 is given by 


'^grav 


2GMb I 2GMb sin sin 0 cos A 

COS^-h ^^^3 


2GM 


‘ grav 


In 


I -h 


-(r 2 • b)/r 2 - (b • s) + ^ - 


ri + (r 2 • s) 


(3) 


■ GM b^ 

c3 r|(l-cos0) 

_ GM b^jcos (f cos ^+sin ip sin 6 cos A)^ 

c3 r2(l-cos0) 


( 8 ) 


Lormula ([^ may be re-written using the approximation ln(I -l- 
y) ^ y -y 1^ iov y « I because the baseline length is much 
smaller than the barycentric distance to the points on Earth 


2GM 

r3 


-(r2-b)/r2-(b.s)+^-'^|^ 

r2+(r2-s) 


1 / -(r2-b)/r2-(b-s) \^ 

2 V r2+(r2-s) ) 


(4) 


One has to keep the barycentric coordinates r 2 for the second 
station in the formulae below to provide consistency with the 
original gravitational delay model 0. 

Converting the three dot products with angles p, ij/, 6 
(Lig[T]), where 

(b • s) = |b| cos^, 

(b-r2) = |b||r2|cosj/r, (5) 

(r2-s) =- |r2|cos^ 


GM b^ sin^ ip sin^ 9 cos^ A 
c3 r\{\-co^9y 


GM b^ cos^ (f GM b^ sin 2ip sin 9 cos A 
c3 r\ c3 r2(l-cos0) 

Using the small-angle approximation {R « r 2 ), 
cos 6^1-^ and sin ^ -, 

2r| r 2 

where R is the impact parameter and is equal to the solar radius 
for the case of grazing light. Therefore 

b%m9 ^ 2^ b^ ^ 2^ b^ sin 9 ^ 2^ /qx 

r2(l-cos0) R ’ r2(l-cos0) ’ r^Cl-cos^) ^2^’ 

and with a 'typical' baseline of ~ 6,000 km the first term ~ 160 
nsec, the second term is ~ 1.2 ns, and the third term is ~ 3 ps 
for grazing light. As we show below, the third term is less than I 
ps for observations at ^ > 1° from the Sun. This term and other 
small terms in ^ may be neglected because VLBI observations 


Article number, page 2 of[^ 












































Titov O. and Girdiuk A.: Applying the theory of general relativity to reducing geodetic VLBI data 


are not undertaken closer than 4° from the Sun. Then, the gravi¬ 
tational delay is given by 


grav 


2GMb __ , 2GMb sin sin 0 cos A 

r2C^ ^ r2C^ 1-COS0 


■ GM b^ _GM b^ cos^cpcos^ 6 

r^il-cosO) r^il-cosO) 

_ GM b^ sin^ (p sin^ 6 cos^ A 
r|(l-cos0)2 ’ 


or, finally. 


'^grav 


2GMb 


cos^ -f 


2GMb sin ip sin 6 cos A 
r2C^ 1-COS0 


, GM ^^(1-cos^ ^cos^ 6) 
r^Cl-cos^ 


GM b^ sin^ p sin^ 6 cos^ A 
r|(l-cos0)2 


( 10 ) 


( 11 ) 


This equation is sufficient for a picosecond level of accuracy. 
The first term of Eq. in) can be excluded because it cancels out a 
similar term from the total geocentric delay ( [Petit et al.|2Q10| ) Eq 
(11.9). Indeed, this formula (11.9) from the lERS Conventions 
includes the term 



Fig. 2. Gravitational delay coordinate term, and main deflection 
term |T3] ) with respect to the angle 6 (see more details in the text). 


Eor an arbitrary 6 the gravitational delay Q can be split into 
several terms in Eq. and the main defiection term is given 
by 


_ (7 + l)GM (b • s) 

coord — ^ 

c^r c 


( 12 ) 


Recalling that y= 1, the resulting effect of general relativity, 
tg/?, in the total VLBI delay model is therefore given by the sum 
of the gravitational delay Tgrav and the coordinate term t coord 


2 GMZ 7 sin ^ sin ^ cos A b . ^ 

tgr = -^- 7 — = 0(2 • - sm^cosA, (15) 

c^r 2 1 - cos^ c 

where a 2 is the defiection angle at arbitrary 6 ( |Shapiro| 19671 
|Ward|197'0| ), and is given by 


Tgr 


r +T , - T + 

‘ grav ^ ‘ coord — ‘ grav ' ^2 ^ ^ 


Tgrav + COS (f. 


(13) 


The gravitational delay ^TT) includes the barycentric dis¬ 
tance of the second station r 2 in the denominator, whereas the 
term includes the barycentric distance of the geocentre r. 
Eor the case of the Sun, the coordinate term may be as large as 
400 ps (see Sect. and should be considered carefully, but 
for Jupiter it does not exceed 1 ps. 

The difference between two terms — (fl^ and — 

_ c^r c Y f c^r2 c 

( pA] ) is less than 0.1 ps, and it may be ignored for me current level 
of accuracy. Nonetheless, the angles 6 and if/ are to be referred to 
the vector r 2 for correct calculation. 

In the original form the gravitational delay is expressed in 
the BCRS system, whereas the baseline length in ( p^ is given 
in the GCRS system. This means that theoretically, one has to 
convert the baseline length h and r 2 into Eq. ( pA) from BCRS to 
GCRS before comparing the coordinate terms from two different 
parts. However, in practice, the corresponding correction in time 
delay does not exceed 1 ps for grazing light and a baseline of 
6,000 km. 

After all reductions, the total contribution of the relativistic 
effects to the total VLBI delay, tgr, can be written as follows: 


2GM sin 6 
c^r2 1 - cos^’ 


(16) 


The minor terms in may be ignored for a sufficiently 
large 6. Therefore may be used as an alternative to the con¬ 
ventional total delay formula for all sources that are far from to 
a gravitational body. Thus, formulae ( p3] ) and ( p^ present the 
defiection angle in a more convenient form for estimation from 
geodetic V LBI data than previously ( jTreuhaft & Lowe|1991t|Ka-| 
|plan|1998] ). These formulae also prove that the two approaches 
- based on angular observations and on the VLBI group delay - 
are equivalent for all ranges of elongation angles, therefore the 
difference between two types of quantities is equal to zero. 


2.1. Alternative equation for the gravitational fields of the Sun 
and Jupiter 

We consider two realistic cases and compare the alternative 

Tgrav Q- There 


formula tgr (14) with the conventional formula 


are two bodies in the solar system that are able to cause a 
significant impact on the observations - the Sun and Jupiter. The 
Sun is the more massive of the two bodies, but observations of 
radio sources closer than 4° are very rare. Jupiter is three orders 
of magnitude less massive, but observations at its limb (i.e., at 
an extremely small impact parameter of 70,00 0 km) are pos - 
sible. Other massive planets have a weaker effect (|Klioner|2003 ). 


Tgr 


2GMZ? sin^sin^cosA GM ^^(1-cos^ ^cos^ 0) 
r 2 C^ 1-cos 0 c^ r^(l-cos6) 


GM b^ sin^ p sin^ 6 cos^ A 
c^ r^il-cosOy 


12 h. 


The alternative formula (H consists of three terms: the 
first term, ti, provides the main contribution to the GR effects, 
^ ^ while the two other terms, t 2 and ^ 3 , contribute much less than 

h. However, we show that for very small angles {6 less than 
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30") all three terms have a similar contribution. At ^ = 180° 
the relativistic contribution, tgr, to the total delay is equal to 
zero, whereas the gravitational delay, Tgrav, is up to 400 ps at a 
'standard' baseline of ~ 6,000 km as a result of the coordinate 


term (12) - 


2GMb 

c^r2 


COS^. 


We illustrate this in Fig. where the coordinate term ( p^ 
and the main deflection term (first term in ( p^ ) are presented as 
a function of the impact parameter for A = 180° (cos A = 1) and 

(p = 45° (cos ^ = sin ^ b = 6,000 km, and r = 1 a.u. The 

coordinate term is always equal to 280 ps because it does not 
depend on the variable impact parameter 0. The gravitational 
delay model Q does not contain the angle 6 explicitly, the solid 
curve in Fig. ^is calculated numerically. It is equal to zero at 6 
= 90° and -280 ps at ^ = 180°. The main deflection term, ti in 
, is equal to zero at ^ = 180, as was previously noted for the 
relativistic contribution, tgr. 



To investigate the accuracy of the alternative formula ( pA) , Fig. 4. Contribution of the minor terms to the total GR delay for differ- 
one has to sum the original gravitational delay formula Q and ent baseline lengths, 
the term ( p^ 


2GM ^ Iril-r(s-ri) 2GMb 
tgr 2 = —^ In r-^—;-T ^ cos ip 

I- I . - X 


(17) 


|r 2 | + (s • r 2 ) 

Figure shows the difference between the delays tgr and 



Fig. 3. Difference between formulae IE) and (It) for a baseline of 
10,000 km length. 

tgr 2 ( p^ for the Sun and a baseline of 10,000 km: the differ¬ 
ence between the two formulae only exceeds 1 ps at ^ = 0°.5. 
Thus, the difference between the two formulae is negligible for 
the light propagation in the solar system as measured by geodetic 
VLBI at the current level of precision. 

While the two minor terms, t 2 and contribute less than 
ti , they cannot be ignored completely. The contribution of these 
terms grows rapidly with baseline length (t - b^) and is 72 ps at 
e=rfovb = 10,000 km (Figg. 

Herewith, we consider two observational VLBI sessions to 
separately investigate the effect of 0 separately during a standard 
24-hour session for the Sun and Jupiter. The VLBI experiment 
R&D1208 was carried out on 2-3 October, 2012 to track sev¬ 
eral reference radio sources in close vicinity to the Sun. While 
the Sun travelled along the ecliptic during the 24-hour session. 
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quasar 1243-072 was tracked in the range of angular distances 
from of 3°.7 to 4°.3. The four curves in the left-hand of plot of 
Fig. reflect the variation in ti for four selected VLBI base¬ 
lines of differing length: Kokee - Tsukub32 (KT; 5,755 km), 
HartRAO - Wettzell (HtW; 7,832 km), OnsalabO - Wettzell 
(OW; 920 km), and HartRAO - OnsalabO (HtO; 8,525 km) with 
respect to the angle 6 of the approach of the Sun to the radio 
source 1243-072. This varies steadily due to the relatively slow 
apparent motion of the Sun. 




arcseconds 


Fig. 5. Main term q for 1243-072 (left) and for 1922-224 (right) vs 6. 


The largest solar system planet Jupiter approached the radio 
source 1922-224 on 18-19 November, 2008. This event was ob¬ 
served during the 24-hour experiment OHIG60. Four stations - 
Hobart26, Kokee, Tsukub32, and Parkes participated in track¬ 
ing the radio source over only 12 hours because of observing 
constraints. The angle 6 between the radio source and Jupiter 
changed from 1'.4 to about 5' within 12 hours. The right plot 
of Fig. shows the variations in q as a function of 6 during 
this event. This term reaches its maximum of about 400 ps for 
the longest available baseline (Hobart26 - Tsukub32 (HT) of b 
=8,088 km) and becomes negligible after a short period of time. 
For the shorter baselines Parkes - Tsukub32 (PT; 7,233 km) and 
Hobart26 - Parkes (HP; 1,089 km)), this term is lower in magni¬ 
tude. 

The coordinate term calculated for both events (Fig.[^ is sig¬ 
nificant for the Sun (reaching 300 ps) but is negligible for Jupiter 
(less than 0.1 ps). In accordance with Eq. ( p^ , the coordinate 
term does not depend on the impact parameter and, therefore, it 
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Fig. 6. Variations of the coordinate term for 1243-072 (left) and for 
1922-224 (right). 




Fig. 7. Comparison of the two models for the approach of the Sun to 
1243-072 (left) and of Jupiter to 1922-224 (right). The difference for 
baseline Hobart26 - Parkes in the right plot is not distinctive for this 
scale. 


does not grow at small angles, as seen in Fig.|^ All variations in 
Fig.j^are caused by the variations in angle ip and baseline length. 

Figure [7] illustrates the residuals between Eqs. ( p^ and ^Tf) 
for these two events. The difference does not exceed 1 ps for ei¬ 
ther case. This proves that Eq. ( p^ , proposed in this paper, can 
be used to model the relativistic effects along with the conven¬ 
tional model ( p^ . 


3. Small-angle approximation. 

The monopole effect of the general relativity results in a circu- 
lar motion of the d eflected object ( [Sazhin et al.|200lt |Kopeikin| 
& Makarov||2007| ). Therefore, the radio sources draw a regular 


pattern on the sky due to the light deflection. Here we consider 
an impact of the minor terms in Eq. ( p^ on the formula for the 
light deflection, in particular on the intraday time-scale, if a large 
planet such as Jupiter passes by a reference radio sources very 
quickly. In the small-angle approximation, Eq. the effect of 
GR ( p^| ) is given by 


Table 1. Einstein deflection angle and the secondary deflection angle 
for the Sun. 


Ri, km 

e 

^2 

^2 ’ 

b = 3,000 km 

mas 

b = 10,000 km 

700,000 

0°.25 

1".75 

3.75 

12.5 

2,800,000 

r.o 

0".44 

0.23 

0.80 

5,600,000 

2°.0 

0".22 

0.06 

0.20 

14,000,000 

5°.0 

0".09 

0.01 

0.03 


Table 2. Einstein deflection angle and the secondary deflection angle 
for Jupiter. 


R 2 , km 

6 

mas 

^2 ’ 

b = 3,000 km 

mas 

b = 10,000 km 

70,000 

20" 

16 

0.34 

1.14 

210,000 

1 ' 

5 

0.04 

0.13 

420,000 

2' 

2 

0.01 

0.03 






Fig. 8. Arc for baselines Hobart26-Tsukub32 (upper left), Kokee- 
Tsukub32 (upper right), Hobart26-Parkes (bottom left), Kokee-Parkes 
(bottom right). 


length in the small-angle approximation as follows: 


^GR = ^ ^ sin ^(cos A - ^ sin ^ cos 2A) = 

= 0 ^ 2 ' ^ sin ^(cos A - ^ sin ^ cos 2A), 


a = — 

Z? sin ^ cos A 

(18) 


4GM /1_ b cos 2A sin ip v 

2/?2 cos A I 

,/ /1 b cos 2A sin (p v 

^2^ “ ^ I 


( 20 ) 


where is the deflection angle for light propagated through a 
gravitation held ( |Einstein|1916| ) for the second station 


a 


2 “ 


4GM 


(19) 


— (x'2 + (^2 • 

Thus, once in single-telescope mode, only the main angle 
is presented; in the two-telescope mode, the secondary de¬ 
flection angle 


and R 2 = r2COS^2- Re-arranging formula ( p^ provides an ex- 
pression for the deflection angle, for an arbitrary baseline ^2 


2GM b cos2Asin^ 
c^R 2 Ri cos A 


( 21 ) 
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represents the difference between angles and a [, if a[ = 
and Ri = n cos 6 i. Even though the parallactic shift of the Sun 
6 i - 62 is about 8" for b = 6,000 km, this difference between a\ 
and 0^2, is not negligible in the small-angle approximation. The 
formula shows that the secondary deflection angle depends on 
the baseline length, but it reduces to Einstein’s classical expres¬ 
sion for a zero-length baseline (i.e. h =0). 

Table H shows the meaning of the angle a'^ for the Sun and 
Table [^fo^upiter with different combinations of the parameters 
O 2 and R 2 . This effect rapidly grows as h is decreasing and may 
be detected with a long-baseline radio-interferometer at a very 
close approach of the Sun (< 5°) or large planets (< 1') to a 
distant radio source. 

The angle a'^ can change comparatively fast when Jupiter 
passing by a radio source. An example of the fast variations of 
a 2 is shown in Eig. [^for the appr oach of Jupiter and the radio 
source 1922-224 discussed in Sect. |2.1[ The corresponding time 
delays as shown in the right plot of Eig. 5 are different, but the 
light deflection angles in Eig 8. for all baselines are the same. 


4, Post-post-Newtonian effect 


As shown by ( |Klioner||l991| ), the main term of the post-post- 
Newtonian effect is given by 


TPPN - 


4G^M^ 


c 5 Llril+hi-s) |r2l+(r2-s) 

Using the mathematical transformation of Eqs. and ([^, 


TPPN - 




r2(l-cos0)| 


b cos ij/+b cos if \ 
^ r2(l-cos0) j 


r2(l-cos0) 


( 22 ) 


(23) 


If the baseline length is negligible with respect to the barycen- 
tric distance (b « r 2 ), one could recall x = ^nd 

apply the approximation = (1 -h v) using Eq. M to solve the 
spherical triangle (Eig. 1) and to develop as fellows: 


TPPN 


4G^M^b -cosip cos 0-sin ip sin 9 cos A+cos ip 


4G^M^b cos(^ 

-cos 6) 


4G^M^b sin ^ sin 0 cos A 
r|(l-cos0)2 


(24) 


In the small-angle approximation {R « r 2 ) (24 (is given by 


^ _ 4G'^M^bcos(p l6G^M^b r2sm(pcosA 

TPPN - —^-R2 - 3 - Ri -• 


(25) 


For the grazing light the first term in Eq. is about 1 ps, 
whereas the seco nd one can rea ch 300 ps at a ’standard’ base¬ 
line of 6,000 km ( |Klioner|1991| ). 

When the first term in Eq. ( [24| ) is negligible, we can intro¬ 
duce a new variable 


_ 4G^M^ 

O^PPN - - 


sing 


(26) 


r^Cl-cosgp 

and re-write the second term of Eq. p4| ) as follows: 

TppN = sin^cosA. (27) 

It is obvious that the variable appj^p plays the role of the post- 
post-Newtonian deflection angle at arbitrary 6 , and converting 
between the deflection angle and the corresponding part of grav¬ 
itational delay uses the same factor - sin <p cos A for the post- 
Newtonian effect as in Eqs. ( p3] ) and ( pfe| ). It enables us to present 
the sum of the two delays tqr ( [15] ) and tpp^^ ( [27] ) in the form of 
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two angles, Eqs. (19) and (26), multiplied by ^ sin ^ cos A, that 
is. 


^GR + '^PPN - ip'GR + O'ppn)- sin^cosA. 

c 


(28) 


Einally, we recall that the second term of Eq. (24) corre¬ 


sponds to formula (11.14) from the lERS Conventions (Petit et 
|al.|2010D . 


5. Proper motion of extragalactic radio sources 
induced by the Sun 


The solar gravitational field will result in a proper motion of ex¬ 
tragalactic radio sources if the effects of general relativity are not 
accounted for in VLBI data reduction. The rapid change of the 
observed positions during a typical 24-hour VLBI session may 
be detected for observations within 10° from the Sun. 

The partial derivative is given by 

drcR 2GMsin ^ cos A 
do r(l-cos^) 

Eor an arbitrary angle 6 , the proper motion p = ^ (where 
a is from ([Tg) may be presented by re-arranging Eq. ( |29| ) as 
follows: 


da c drcR c drcR dO 

p = — ~-=- 

dt Z? sin ^ cos A dt Z? sin ^ cos A dO dt 

2GM 1 do 
c^r (I - cos 0 )dt 


(30) 


Eigure [^illustrates Eq. ( [30j ) for a range of angles 0 from 0° to 
10°. It is obvious that all distant radio sources across the sky will 
display a circular motion if general relativity is not accounted 
for in VLBI data reduction. The magnitude of the effect is pro¬ 
portional to . In particular, a source near the ecliptic poles 
0 = 90° will draw a circle of 4 mas in size. The daily displace¬ 
ment between the consecutive days will be equal to ^. Eor 


the Sun, ^ = 0.0174, and the daily displacement is 0.07 mas. 


Eor all other sources, the daily displacement will be proportional 


to the factor 


1 


(l-cos 6) 


in Eq. (30). 



Fig. 9. Daily proper motion calculated from Eq. { [^ near the Sun. 


In the small-angle approxiation, Eq. @ may be written as 
follows: 
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4GMr de _ r de 
~dt ~ ~^R~dt' 


(31) 


For background radio sources that are co-planar with the 
plane of the body, the apparent motion, yu, rapidly grows from 
0.5 yuas/year at ^ = 10°, to 1.6 yuas/year at ^ = 5° and 40 yuas/year 
at ^ = 1°. Therefore, a systematic pattern of proper motion in a 
particular area of the sky may serve as an indicator of a hidden 
mass in this direction. 


6. Conclusion 

The effects of general relativity are explicitly contained in both 
components of the total VLBI delay model - the gravitational 
delay and the coordinate term with a factor of 1 - . While the 

former component uses the individual barycentre positions of the 
radio telescopes to calculate the effect, the latter component is 
expressed in terms of the baseline between the radio telescopes. 
Coupling between both parts has not been investigated until now. 

The Shapiro effect as measured by radars and the deflec¬ 
tion of light measured with traditional astronomical instruments 
are considered two independent tests of general relativity. We 
showed that the total group delay model joins the two tests 
within one observational technique - geodetic VLBI. The grav¬ 
itational delay that traditionally originated from the Shapiro ef¬ 
fect is linked to the light deflection angle. Therefore, the two 
approaches, VLBI delay and angular, are absolutely equivalent. 

For almost all realistic situations this angle does not depend 
on the baseline length, thus, a standard geodetic VLBI interfer¬ 
ometer acts as a traditional astronomical instrument. In addition, 
the coordinate term explicitly presented in the conventional to¬ 
tal delay model ceases to exist because it is compensated by the 
same effect in the gravitational delay with opposite sign. Thus, 
the proposed alternative version of the general relativity contri¬ 
bution to the total VLBI group delay model is free of the co¬ 
ordinate effects. Therefore, the two approaches, time delay and 
angular, are absolutely equivalent. 

The final equation of the general relativity contribution ex¬ 
pressed in angular terms also comprises two smaller terms that 
are significant at very small angular separation between the de¬ 
flecting body and distant radio source. These terms may be con¬ 
sidered as an increment in the light deflection angle due to the 
additional time delay for the propagation of the light from station 
1 to station 2. This effect becomes significant at | >0.1. 

It is unlikely that the effect will be detected with ground- 
based interferometers, but it may be detected with space-baseline 
interferometers, for example. RadioAstron with b up to 300,000 
km ( [Kardashev et al.|20T3] ). 
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